General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



. ,r 



A IiX:AL EARTHQUAKE CODA MAGNITUDE AND ITS RELATION 
TO DURATION, MCMENT AND LOCAL RICHTER MAGNITUDE Mj^ 


By 


Anne M. Suteau 
James H. Whitcomb 


Seinmological Laboratory 
California Institute of Technology 
Pasadena, California 91125 


(NASA-CR-157473) A ICCAL EASTHQUAKE CODA 
HAGNITUDE AND ITS RELATION TO DURATION, 

HOHENT M SUB O, AND LOCAL RICHTER MAGNITUDE 
H SUB L (California Inst, of Tech.) 35 p HC 
A03/MF A01 CSCL 08K G|446 


N78-29664 


Dnclas 

27244 


AUO 1976 


UVi 






Vi?, 




■\y 



ABSTRACT 


A relationship was found between the seismic moment, of shallow 
local earthquakes and the total duration of the signal, t, in seconds, 
measured from the earthquake origin time, assuming that the end of the 
coda is composed of backscattering surface waves due to lateral hetero- 
geneity in the shallow crust following Aki. Using the linear relationship 
between the logarithm of and the local Richter magnitude M^, we 
obtain a relationship between and t, of the form: ^ i 

1/3 

a 2 t + f(t), where a^, a^^, 32 constants depending on Q and geometric 

spreading, and f(t) is a function of the instrument response and a (weak) 

function of the scattering process. This relationship is different from the 

2 

(log t) + a^A, where T 
is the duration measured from the first P arrival time and A is epicentral 
distance in km. In the theoretical relationship, the dependence on epicentral 
distance is implicit in t. This relationship has been used to calculate 
a coda magnitude which has been compared to for Southern California 
earthquakes which occurred during the period from 1972 to 1975. This 
comparison has been made independently at six stations of the CIT network. 

M^) has been obtained. 

The standard errors range from .2 to .3 and the correlation coefficients 
from .80 to .90. Once station gain is accounted for, station correction 

and 

estimation is not limited to a duration measurement but can utilize the 
entire earthquake coda. The use of several measurements per seismogram 
produces a more robust estimate of earthquake size. 


terms are less than 0.17 magnitude unit when comparing 


At all stations, a good linear fit (M^^ ^ ^ ^1 


empirical one generally used ^ ^1 ^ ^2 
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INTRODUCTION 

With the advent of computerized local earthquake analysis systems 
using large seismographic networks, new methods of earthquake magnitude 
determinations are needed. Current methods of magnitude determination 
for local earthquakes, estimations that require the amplitude of body 
waves or the duration of the coda, are not generally useful for computerized 
systems. In larger earthquakes, body wave amplitudes can exceed the linear 
range of the recording system and, in some automated systems, recording 
can cease before the end of the coda. The intent of this paper is to 
develop a method of magnitude estimation that can utilize any part of 
the coda, some of which will certainly be preserved in the linear range 
of computer recording systems. It is shown how this coda magnitude method 
relates to other estimates of earthquake size by the use of coda duration, 
moment and local Richter magnitude 

An empirical correlation between the magnitude of earthquakes, M, and 
the duration of the recorded signal, t, has consistently been observed. 
Linear relationships between M, log t (t in seconds), the epicentral 
distance, A (in km), and the focal depth, h (in km), have been empirically 
determined for various geographical areas, using multiple regression analy- 
sis techniques. 

The first study was made by Bisztricsany (1958), who found a linear 
relationship between the magnitude of teleseisms (in the range of 5 to 8) , 
the logarithm of the duration of the surface wave train, and the epicentral 
distance. Later on, this method has been applied to local earthquakes 
(of epicentral distances less than 100 km), with shallow focus (focal 
depth less than 60 km), and the duration T has been defined as the total 
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length of the record, instead of the length of the surface wave train 
only. In chronological order, such studies have been due to Solov'ev 
(1965) for Sakhalin Island; Tsuraura (1967) for Japanese earthquakes 
recorded at the network of the Wakayama Microearthquake Observatory, 
and whose magnitudes were determined by the Japan Meterologlcal Agency; 
Crosson (1972) for earthquakes In the Puget Sound Region, using magni- 
tudes calculated from Wood-Anderson records, Lee et , (1972) for 
central California earthquakes with Richter magnitude M|^; Real and 
Teng (1973) for Southern California also with M^; Hermann (1975) for the 
central United States and body wave magnitude determinations, Bakun and 
Llndh (1976) for the Oroville, California, region, and Similar 

studies have been made for the Mica Array, Canada, the Santiago station 
in Chile, the Alaska network and the University of Utah seismograph 
stations, U.S.A. , etc. ...(see Lee and Wetmlller, 1976). 

For shallow earthquakes, focal depths are much smaller than epicentral 
distances and much less accurately determined. As a result. It Is harder 
to Investigate the dependence of duration on focal depth than on epicentral 
distance. No such dependence has ever been found in a conclusive way. 

The local magnitude M^^, Introduced by Richter (1935), Is given by: 

M^ > log A - log A^ (A, km) where A Is the maximum amplitude In mm on the 
Wood-Anderson seismogram for an earthquake at distance A km, and A^ Is 
that for a particular earthquake selected as standard. (All logarithms used 
here are common logarithms to the base 10.) The magnitudes determined 
by this method are usually above 2.5, so that, for smaller earthquakes, 
an equivalent magnitude has to be used. It Is obtained by reading the 
maximum amplitude recorded by a short-period high gain seismometer and 
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conversion into the amplitude that would have been recorded if a 
Wood-Anderson Instrument had been used. 

When plotting the local magnitude versus the logarithm of duration » 
a slight curvature Is consistently observed (Lee , 1972, Figure 4, 

Real and Teng, 1973, Figures 5 and 6; Bakun and Llndh, 1976, Figure 4a; 

Lee and Wetmiller, 1976, page 23). An example for the station MWC is 
shown in Figure 1 from Real and Teng (1973). The slope increases with 
increasing magnitude. Consequently, a better fit of the data is obtained 
by introducing a quadratic term in log t in the magnitude — duration 
relation. Therefore, the general form of the empirical relation between 
magnitude and duration is given by: 

M » a^ 4- a^^ log T + a 2 (log t)^ + a^A, (1) 

where a^, a^^, a 2 and a^ are constants. Reviews of existing formulas 
have been given by Real and Teng (1973) and Bakun and Lindh (1976). The 
standard deviations are in the range of .1 to .2, the estimated 
standard errors, of .2 .3. The dependence on epicentral distance is 

small and always negligible for distances less than 200 km (e.g., Tsumura, 
1967; Crosson, 1972). 

As established by Tsumura (1967) , the errors in the duration magni- 
tude due to subjective errors in reading durations do not exceed ± .3 magni- 
tude unit. The main source of observational error is the uncertainty in 
determining the end of the earthquake signal. Accw*'ding to Real and Teng 
(1973), it is fixed with an uncertainty of about 20%. It results in an 
approximately constant uncertainty of log t over the range of magnitude 
of interest, .15 unit of magnitude. In the same study, it has been 
established that there is no azimuthal dependence, as long as a sufficient 
number of events covering a large area and recorded on a long period of 
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time (say, a few years) Is analyzed, and that the effect of variations 
In mlcroselsmlc noise is Insignificant in Southern California. 

The >retical Relationship Between and Duratio n 

In this paper, we try to find a theoretical Justification for the 
empirical relationships existing between local magnitude, duration of 
the signal, and eplcentral distance. Akl (1969) has proposed a theory 
in which the coda of local earthquakes is composed of backscattered waves 
due to lateral heterogeneity in the crust. He has observed that the 
power spectra of coda waves at a given time measured from the earthquake 
origin time ("lapse time") appear to be nearly Independent of the 
eplcentral distance and of the nature of direct wave path between station 
and epicenter, as can be seen by comparing the records of the same event 
at different stations (Akl, 1969; Aki and Chouet, 1975). Aki (1969) 
has constructed his model for coda waves on the following three assumptions: 
(1) The distribution of the scatterers is two-dimensional over the earth's 
surface, random, and uniform. (2) The primary waves and the secondary 
waves are surface waves of the same kind, and a constant group velocity 
Independent of frequency is assimed for them, for the sake of simplicity. 

(3) The distance between station and scatterer is on the same order as 
the distance between hypocenter and scatterer and is much larger than the 
hypocentral distance. He makes use of the Born approximation by neglect- 
ing multiple scattering. He obtains the following relation: 

M^[2N(r^)]^''^U^(fp|r^)| - t^^^ 


( 2 ) 
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where M is the seismic moment, N(r ) Is the number of scatterers within 
o o 

a radius of r , U (f |r )| Is the absolute value of the Fourier Transform 
o ' o p' o ' 

of displacement due to secondary waves generated at a scatterer at a dis- 
tance r_ by a source of unit moment located at the same distance from the 
o 

scatterer, A(t) Is the average peak-to-peak amplitude on a seismogram estimated 

from a short time sample around a given absolute lapse time t, I(fp) is the 
Instrument response, and Q Is an apparent quality factor, which is assumed 
constant and Includes the effects of Intrinsic attenuation of the wave 
medium and scattering process. The peak frequency f^ was obtained by 
measuring the wave period directly on the seismogram at time t. For 
his data (aftershocks of the Parkfield earthquake, June 28, 1966), 

Aki found the following empirical formula for the peak frequency f^ versus 
the corresponding lapse time t: 



where t is In seconds and f is in cps, in the range 30 < t < 1000. This 

P 

relation seems to be Independent of earthquake size, eplcentral distance, 

and nature of the direct path between station and epicenter. In the 

following, we will assume that the same dispersion relation applies to 

earthquakes studied here. Let B(f ) = [2N(r (f |r )|. The 

p o o p o 

values of B(fp) as a function of peak frequency have been empirically 
determined by Aki and are given in his Table 3 (Akl, 1969). Using (3), 
(2) becomes: 


1 / 3 , 


M » B(f )“^ ^11/12^(1^17.67 t /Q) 1 14 


5.24 


licy 


(4) 


• A(t) 


-1 
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This formula is strictly valid only in the range of frequencies which 
are lower than the corner frequency of the particular event under consid- 
eration. The effect of the radiation pattern has been neglected, which 
is a reasonable assumption because the scattered energy should represent 
a broad range of source azimuths. A linear relationship between the 
logarithm of seismic moment and the local magnitude for a given 
geographical area has been consistently established in California over a 
wide range of magnitudes (0 < < 6) (e.g., Aki, 1969; Wyss and Brune, 

1968; Thatcher and Hanks, 1973; Johnson and McEvilly, 1974; Bakun and 
Lindh, 1976). No systematic dependence of this relationship on 

either epicenter location or focal depth has ever been found* Uncertain- 
ties in the determination of magnitudes and variation In the geology at 
the recording site introduce an uncertainty factor of 2 or 3 in relating 
moment to magnitude (Wyss and Brune, 1968; Johnson and McEvilly, 1974). 

We will express this relationship in the following way: 

log = a -h b (5) 

where a and b are constants. If we designate as the coda magnitude 
Mq, (4) can then be rewritten as: 

a + bM- * -0.72+ log A(t) + log B(f ) ^ 

^ P (6) 

1/3 

- log Q + ^ log t + ^ 17.67 log + log |l(fp) 

Then we divide the term corresponding to the instrinaent effect Into 
two parts: the gain G, and the normalized amplitude 


spectrum R(^p) • We now have the relation for coda magnitude 


"c ■ 


a - 0.72 + log^Vlog R(f )~^ 

G p 


+ log B(fp)"^ - log Q + -jl log t + t 


1 / 3 ^ 


(7) 
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Two different definitions of signal duration x are coinnonly used. The 
time is always measured from the onset of the first P arrival. The 
difference occurs in the way of defining the end of the coda. It has 
been defined as the time when the signal returns to the noise level 
present before the arrival of the P wave (e.g.« Crosson, 1972) (Criterion 1) 
or as the time \idien the trace amplitude (peak-to-peak) falls below a 
given level on some reading device (e.g.« Lee e^ , 1972) (Criterion 2). 

If t in expression (7) is the time when the coda ends, measured from the 
earthquake origin time, we have: 

t - + T (8) 

a 

where a is nearly a constant depending on the P wave velocity* With 
Criterion 2, A in expression (7) is a constant. With Criterion 1, A 
represents the average amplitude of the background noise* It is a func— 
tlon of the recording station, but the ratio - is nearly constant for all 
stations in a given area, since an increase in the magnification of the 
instrument results in an increase in the average noise level* 

Hermann (1975) has studied the Influence of different cut-off levels 
for defining duration for earthquakes with magnitudes between 2*5 and 4*5 
on the magnitude-duration relation. Duration measurements made at two 
different levels at each station yielded relationships differing from 


one another by only an additive constant and, within the limits of accuracy, 
the linear relation between magnitude and log t were equivalent at the two 
levels, as expected* In the following, we will consider A as a constant* 


(7) can be rewritten in the following way, when using (8); 


- r F- a - 0-72 + log ^ + log R(f )“^ + log B(f ) 

C b L ^ P P 


-1 


- i log Q + II log (| + t) + (^ + T)^^^]. 


( 9 ) 
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In agreement with the empirical relationship (1) between t, and A, 
is found to be an increasing function of both t and A, but 
in (9) is more complicated. In addition to the logarithmic term in 
duration, there is a power term; and the effect of the gain of the instru- 
n^nt has already been mentioned. For a given area, Q has been considered 
as constant. Aki (1969) fixed it at 200. The term log R(fp) ^ depends 
on the instrument type. As it is an implicit function of i and A 
(through (3) and (8)), the instrument is likely to influence all the 
coefficients in the empirical formula (1). This result is consistent 
with the observations of Bakun and Lindh (1976): They plotted vs. 

log T for three different instruments and found a different Intercept 
and a different slope in each case (their Figure 4b). The USGS and ORV 
instruments, which have about the same normalized response below 5 hz, gave 
about the same slope, whereas the ORV-LP instrument which has a normalized 
response very different from the two others gave a different slope (their 
Figure 2) . 
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DATA ANALYSIS 

In order to further check the duration-M^ relationship, we now 

compare local magnitudes with from the theoretical expression (9) 

for events in Southern California recorded during the period from 1972 

to 1975. The following stations of the telemetered seismic network 

operated by the Seismologlcal Laboratory of the California Institute of 

Technology have been used: Isabella (ISA), Goldstone (CSC), Saddleback 

Butte (SBB), Mount Wilson (MWC) , Twentynine Palms (TPC) , and Palomar (PLM) . 

Their locations in the Southern California Seismograph Network are shown 

in Figure 2. The signal was recorded by short-period vertical Benioff 

seismometers and registered on drum recorders revolving at a speed of 60 mm/min. 

Durations t were read using Criterion 2, the average peak-to-peak amplitude 

threshold A being approximately 0.1 cm. The gain, G, was initially assumed 

to be the one at the station PAS, G , of 10^* The effect of different 

gains is discussed below. The coefficient a of expression (7) was taken 

as .162 sec/kro, Q fixed at 200 following Aki (1969) and the values of 

B(f ) were taken from Aki*s Table 3, once f had been determined from x 
P P 

and A using relations (3) and (7). The constants a and b in the moment- 
magnitude relation (5) that we used were also determined by Aki (1969) 
in his: study of aftershocks of the Parkfield earthquake, for the magni- 
tude range of 3 to 5.5 he found: 

log = 15.8 -I- 1.5 M^. (10) 

In this way, we have calculated a coda magnitude for each event at 
each of the si:: stations, as expressed by: 

«c ■ ♦ Sas * 13 

+ ~ log B(f )"^ + 0.61 log (oA+t) + 0.08(aA+t)’^^ 
i . D P 


( 11 ) 
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and shown In Table 1 (for t < 60 sec, values of B(f^) used in Table 1 
are different from Aki, as discussed below). normalized 

Instrument ~ esponse of the short period Benioff seismometer. Fi{;ure 3 
shows the correlation between and ranged from 2.6 to 5.5. 

Below 2.6, most values of have been obtained from short period 
seismograms instead of Wood- Anderson records and they are omitted. 

The results of a least-square fit of the data for each station 
for > 2.6 are given in Table 2, including also the geographical 
coordinates of the stations, the period of time during which the data 
have been obtained, the number of events used in the analysis, the standard 
error of estimate and the correlation coefficient. The standard errors 
are .2 'v .3 which is comparable to what is usually obtained by fitting 
such data with an empirical relationship of the form of expression (1). 

The correlation coefficients range from .80 for MWC where there is much 
scatter in the data, particularly at low magnitudes, to .90 for ISA. The 
number of events is on the order of 300 'V 400, depending on the station, 
and the data have been obtained between 1973 and 1975 for most stations, 
so that the azimuthal coverage is good and the events have a broc.d 
distribution in size, epicentral distance and source mechanism. 


i 
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DISCUSSION 

The relationship between and is put in a ~re convenient form 
when it is written as: 

«L • ♦ <=1 «C 

and these paraiMters are shown in Table 3. It can be seen that the 
slope Cj is systematically too hig^ taking values between 1.22 and 1.55, 
the average value being 1.4. We now evaluate what parameters of the 
theory can affect the slope. 

Aki and Chouet have also proposed a second model for coda waves in 

which the seismic energy transfer was considered as a diffusion process. 

1/2 

With the diffusion theory, the geometrical spreading term t in 

3/4 

equation (2) is replaced by t , i.e., a power intermediate to 1 (corre- 
spending to body waves) and 1/2 (surface waves) . It is easily seen that 
this would change only the coefficient of log (a4 d- t) in the expression 
giving the coda magnitude (11), from .61 to .78. The diffusion theory 
applies only to the parts of the energy scattered in the backward direc- 
tion. Like the single-scattering theory that we have used, it is a 
simplification, and probably, the actual power of the geometrical spread- 
ing term lies between the values of 1/2 (single-scattered surface waves) 
and 3/4 (diffusion process) and the calculated slopes are shown in Figure 4. 
The effect of incorporating a diffusion process theory would reduce the 
observed slopea from 1.4 to about 1.3, which is not a significant difference. 
The coda magnitude for each event recorded at each station has been 
calculated from expression (11). This expression makes use of the empiri- 
cal dispersion law (3) for coda waves found by Aki for the Farkfield region 
(Aki, 1969). A priori, this law may not apply to our set of data and a 
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more rigorous treatment would necessitate a derivation of a similar law 

for the data used in this study. Any discrepancy betii^en the actual 

dispersion law and the one we used may yield a slope different from 1 in 

the vs. curve and even to a curvature of this line. In order to 

get we also used the values of B(f^) tabulated by Aki (1969). These 

values were obtained by fitting the data in his Figure 3. There are not 

many data for peak frequencies between 0.2 and 1 cps, which is the range 

of frequency corresponding to the durations of most events with magnitudes 

between 2.6 and 5.5. Therefore, the function B(f ) is not well constrained 

P 

and may be subject to reevaluation. However, the two most likely candidates 
to change the slope of the vs. correlations are either the 
vs. relationship or Q. 

Relationship 

Wyss and Brune (1968) have Fourier analyzed surface waves from 13 
earthquakes in the Parkfield region and have found the following relation- 
ship between seismic moment and Richter magnitude 


log = 1.4 -f 17.0, 

for 3.2 < < 5.5. of Southern California earthquakes calculated by 

Thatcher and Hanks (1973) from long-period S wave spectral levels yielded 
the relation: 


log = 1.5 + 16.0, 

for between 2.5 and 6. These two studies give a coefficient b of 
1.4 1.5, similar to the one found by Aki (equation 10). However, 

more recent studies in Central California have yielded a coefficient 
b closer to 1. Johnson and McEvilly (1974) have estimated seismic moments 



from the low-frequency levels of tlie spectra for 13 earthquakes with 


magnitudes between 2.4 and 3.1 located near the San Andreas fault and 
found the following relation: 

log - (1.16 t .06)Mj^ + (17.60 > .28). 

And for earthquakes near Orovllle, Bakun and Llndh (1976) have found: 
log - (1.21 1 <17-02 ♦ .07) 

for 0 < < 6. The value of 1.5 for the coefficient of In the moment- 

magnitude relation may be more appropriate for Southern California than 
the value of 1.2 found for Central California. However, a coefficient b 
closer to 1 would exp >in the average slope C^ higher than 1 found for 
our vs. curves. In fact if the actual value of b was 1.2 instead 
of 1.5, this slope would would be reduced from 1.4 to 1.1. 

Effective Q 

In our calculation of coda magnitude M„, we have used a Q 
of 200 (Aki, 1969). However, a lower value may be more realistic, since 
values of 50-200 at frequencies on the order of 1 Hz have been found for 
the Q of coda waves in Central California and Western Japan (Aki and 
Chouet, 1975). At these frequencies (.2 Hz - 2 Hz) which are dominant 
at the end of the coda of local earthquakes, the coda is probably made 
of backscatterlng surface waves from heterogeneities in the shallow, low 


Q crust. 


If the actual average apparent Q of coda waves of local 


earthquakes in Southern California was 70, a value favored by Aki and 
Chouet (1975) t instead of 200, the observed scope of vs. relation- 
ship (12) would be reduced from 1.4 to about 1.1 as seen in Figure 5. 


In order to Investigate the influence of gain of the instrument 
and/or geology at the station site on the coda magnitude 


we now 
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consider the residuals AM S M, 


at the various stations for a 


well-sampled interval of magnitude, chosen around • 3.1, given in 
Table 3. If Is the gain at a station STA, by comparing (9) and 


(11) we get! 


„ . I , STA . „ 

”c ■ r.5 cTT * 


The values of the gain corrections (AG log also 

given in Table 3, as well as the values of AM AG, It is clear from 
this table that the major contribution to the discrepancy between 
and is due to the effect of the gain at each particular station. 

A coda excitation factor, depending on the surface geology of the station 
site, (Akl, 1969) should have the same effect on the magnitude residuals 
as the gain of the instrument. This factor can be 5 to 8 times larger 
on sediments than on granite (Akl, 1969), which corresponds to a rela- 
tive change of ,5 to ,6 In overestimated at ISA by 0,17 

and underestimated at SBB by 0,13, with ttie ottier stations intermediate 
to these values. These corrections may be insignificant considering the 


scatter of the data. All of the stations used in this study are sited 


directly on crystalline bedrock (granite or metamorphic) , with the excep- 
tion of Palomar, which has a few meters of alluvium between the Instruim^nt 
and the granitic bedrock. All are essentially on the ground surface 
except for Isabella, which is about 70 meters in a tunnel in bedrock. ISA 
has an anomalously Irigh residual, AM-AG, of ,17 magnitude unit (Table 3) 
which may be due to its particular setting In the Southern Sierra mountains. 
No t e that the e xc i t a t i on f a c t o r shou I d mo re ge no rail y l^e a f unc 1 1 on o f f re - 
quency , so that it could not only shift the vs. curve but also distort 


it. 


We now compare the theoretical magnitude from equation (9) 

with the empirical curve of our Figure 1 found by Real and Teng (1973) 

2 

for durations at station MWC: « 1.37 + 0.41 (log t) (the distance 

term Is neglected) . If vie assume the same parameters as used above 

in calculation of except for an allowance of variance of Q for a 

good fit, we get the theoretical curve shovm in Figure 1, which is 

calculated with a Q of 500. (Because preliminary results have shown 

no significant variation in the scattering property term B(f^), it is 

-25 

assigned a constant value of 4.32 x 10 which is of course dependent 

on Q. This has only minor effects on Che shape of the low-magnitude end 

of the curve.) In the empirical relationship of Real and Teng, the 

2 

curvature cooes from a (log t) term whereas in the curve it comes 

1/3 a 

from a t term. It is interesting to note that expansion of t in 

a series is: 


“ 1 + ^ t) ^ + ... 

2 

Thus, the next term after log t is a (log t) term in parallel with Real 
and Teng's empirical fit. There is however a discrepancy between Che 
duration data of Real and Teng and that analyzed here. It is represented 
by the fact Chat the former is better fit with a high Q of about 500 
whereas the latter requires a lower Q of about 70. At the present this 
difference is not undet stood. The data of Real and Teng were gathered 
from 1969 to 1971 events and the data here are from 1972 through 1975. 
This presents the possibility of a change of Q between the intervals. 
Alternatively, although Che reading methods were supposedly the same, 
it is possible that variances in definition of Che "end" of Che coda for 
large and small events can occur between readers. Whether these or other 


causes are responsible Is a subject that deserves further study. 

Because the magnitude is a function of instrument response, it 

is interesting to compare the theoretical with duration data of 

Bakun and Llndh (1977) . They found a different duration magnitude 

relation for the Instruments ORV, a Benloff seismometer, and ORV~LP, 

a long period instrument with peak gain at about 0.2 Hz. Figure 6 

shows calculated with a Q of 70 compared with the empirical piece- 

wise linear fits of Bakun and Lindh to ORV and ORV-LP data. (A constant 

-23 

value of B(f ) « 6.24 x 10 was used for the M_ calculation.) The 
p C 

Benloff ORV data are fit well except at the break in the two linear 
empirical curves. But more importantly, the ORV-LP data show markedly 
shorter durations than ORV below * 5.5 in agreement with the theore- 
tical curves. The theoretical curve is not defined below about 

■ 4.4. This implies that, because of the drop off of ORV-LP instru- 
ment response with higher frequency above 6 Hz combined with smaller 
earthquake amplitudes, scattered energy drops below the noise level. 

The development of the theory leading to the expression for in 
equation (9) was based on the assumption that the peak frequency f^ was 
below the characteristic or "corner" frequency of the seismic source. 

This assumption begins to be compromised at values significantly below 
2 because here the characteristic frequency is about 2 Hz (see, for example, 
Johnson and McEvilly, 1974). Further, since the evidence here indicates 
that the duration magnitude scale Is a subset of the general coda magni- 
tude Mg in equation (7) , it follows that extrapolation of duration 
magnitude scales to levels significantly below * 2 is also suspect. 

A further consideration that compromises the extension of duration 
and coda magnitude scales to small earthquakes is the evidence that 
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the later part of the coda Is mainly composed of backscattered waves* 
Certainly this is not true at the time of the direct S-wave arrival* 
Depending on the length of what one considers **S-wave** on the seismo- 
gram, it is reasonable to suppose that within 5 to 10 seconds of the 
S-wave arrival we are dealing more with direct S-wave energy. For 
both this, and the characteristic frequency problem discussed above, 
extrapolations of duration and coda magnitude scales below about » 1*5 
should be carefully reevaluated* 

In Table 1 are summarized the values calculated for the coda 
at t, the time from origin tic^, for 0,1 cm RMS amplitude at a station 
with 1 X 10^ gain. values for Q*s of 500, used in Figure 1, 200, used 
in Figure 3, and 70, used in Figure 6* are given to illustrate the effect 
of variable Q. In these tables, the scattering property B(f^) is assumed 
constant and adjusted such that M values at t = 60 sec are equal. 

Again, the assumption that B(fp) is constant Instead of variable as in 
Aki (1969) significantly affects only magnitudes below « 2.6 and would 
result in a reduction of the lower values by about 0.3 magnitude units 
if the Aki values were used. Study of this parameter will be done in 
a future paper. An Important point to note here is that Q can be an 
Important factor in relating the coda magnitudes of small earthquakes 
to large earthquakes. The difference in coda magnitude can be as large as 
one unit over an range of 3 to 5 for published data in the literature. 
This resalt makes the calibration of individual stations for Q an Important 
consideration in the establishment of a coda magnitude scale. 

The use of values in Table 1 for any part of the coda Involves 
only an additive term for the coda amplitude at the corresponding time. 
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as can be. seen in equation (7). Different Instrument gain requires an 
added constant In the same equation. Use of a different Instrument 
response requires the removal of the Benloff response term In Table 1 
or equation (7), log R(fp) and addition of the equivalent term 
corresponding to the desired Instrument response. 
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CONCLUSIONS 

Backscattering surface waves due to lateral heterogeneity In the 
shallow crust appear to be a good explanation for the late part of the 
coda of local earthquakes and duration times. The scattering 
theory measures the moment and provides a physical basis for the 
observed empirical correlation between local Richter magnitude and 
total duration of the coda. 

With the use of the scattering theory, we derive a coda magnitude 
that can be used at any part of the coda and, in particular, where 
durations are measured. We compare to the coda magnitude calculated 
from duration data at six stations of the CIT network. For a linear fit 
» C^ + Cj^ M^, we find slopes Cj^ systematically higher than one with 
the use of our initial assumptions based on the parameters of Akl 
(1969). This high value of Cj^ Implies one, or more likely a combination, 
of the following factors: the coda waves may have a geometrical spreading 

factor intermediate between surface waves and body waves; the coefficient 
b of 1.5 in the relationship may be too high; the apparent 

Q of 200 that is used is too high and a value of 70 may be more 
appropriate. 

For the six bedrock stations analyzed here, a correction for station 
gain within the theoretical framework of removes all significant site 
correction factors from the duration data. The remaining site correction 
on "coda excitation" terms deviate less than 0.17 magnitude units from 
the average . 

Comparison of the coda magnitude with the duration data of Real 
and Teng (1973) shows that the theoretical matches the curvature of 
the data that led Real and Teng to add a higher order term to their 
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fit of versus log t. However, a high Q of 500 In the relation 

was necessary to fit their data. Cooparlson of with the duration 
data of Bakun and Llndh (1977) shows the effect of different Instrument 
response on the duration data. With a Q of 70, the calculation fits 
their short-period Benioff data from ORV and simultaneously fits the 
widely divergent long period ORV-LP data for magnitudes above “ 4.4, 

The results here strongly support the theory that the coda is composed 
of scattered energy. However, that same theory implies that scattering 
assximptlons do not hold in the early part of the S-wave coda. For this 
reason, it appears that the extrapolation of duration and coda magnitude 
scales below « 1.5 should be carefully reevaluated. 

The coda magnitude M_ can be used either with coda duration data 
or with any portion of the coda that is not "close" to the direct S-wave 
arrival. The latter feature is important for the usefulness of in 
computerized seismic data systems where parts of the coda may exceed 
the system's linear range, preventing the calculation of a body wave 
magnitude, or may be cut off before the time where duration would 
normally be measured. 

Most importantly, the coda magnitude has a physical basis 
that relates moment M^, Richter magnitude M^, duration, and the 
general shape of the coda to the station gain, instrument response, and 
effective Q. There is evidence here that Q can vary significantly from 
station to station which may require a calibration of Individual stations. 
The study of regional variation of Q and other parameters that affect 
the theory of scattered coda energy will be the subject of future 


investigations. 
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TABLE 1 


CALCULATION OF THE CODA MAGNITUDE 
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TABLE 3 


FIT OF + Cj^ 


Station 

* 

Gain 

AG 

C 

o 

^1 

AM 

AM-AG 

GSC 

4.80 E5 

.46 

-1.48 

1.34 

.34 

- .12 

ISA 

2.40 

.26 

-1.73 

1.37 

.43 

+ .17 

MHC 

2.37 

.25 

-1.58 

1.38 

.27 

+ .02 

PLM 

4.83 

.46 

-2.45 

1.55 

.46 

+ .00 

SBB 

8.00 

.61 

-2.10 

1.46 

.48 

- .13 

TPC 

3.76 

.39 

-1.17 

1.22 

.40 
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Figure 2. 


Figure 3. 


Figure 4. 


Figure 5- 


Figure 6. 


FIGURES 

Local Richter magnitude versus log t (duration) for 
station MWC. Dashed line is empirical fit of Real and Teng 
(1973) to their data. Continuous line is the calculation 
of using parameters described in the text. A small 
vertical offset is introduced to separate the curves. 

Locations of the stations GSC, ISA, MWC, PLM, SBB, and 
TPC within the Southern California Seismographic Network. 

Correlation between and for the stations GSC, ISA, MWC, 
PLM, SBB, and TPC. The parameters for the linear least 
square fits are given in Table 2. 

The change in slope of due to different geometrical 

spreading terms (following Aki and Chouet, 1977) for 

1 3/4 

(a) T , backscatter ing body waves, (b) t , diffusion 

1/2 

process, and (c) t , backscattering surface waves. 

The change in slope of due to different values of Q. 

A value close to 70 fits the data of Figures 3 and 6. A 
Q of 5v)0 fits the data of Figure 1. 

Local Richter magnitude versus log t (duration) for 
stations ORV and ORV-LP. Dashed line and squares are 
linear fits to data and data for ORV and ORV-LP, respectively, 
from Bakun and Lindh (1977). Solid lines are the theoretical 
coda magnitude calculated with Q = 70 for ORV and ORV-LP 
to illustrate the effect of different instrument response. 

M^ for ORV-LP is not defined below M^ ~ 4.4 for this 
particular value of Q. 
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